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^ ! Abstract 
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We study pion-nucleon scattering at tree level with a chiral lagrangian of 
' pions, nucleons, and A-isobars using a K-matrix unitarization procedure. 

o\ ' 

Evaluating the scattering amplitude to order Q 2 , where Q is a generic small 
' momentum scale, we obtain a good fit to the experimental phase shifts for 

vo ' 

Q\ ' pion center-of-mass kinetic energies up to 50MeV. The fit can be extended 

Q_i' to 150 MeV when we include the order-Q contributions. Our results are 

independent of the off-shell A parameter. 

43 ' 
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PACS number(s): 11.80.-m, 12.39.Fe, 13.75.Gx 

Pion-nucleon (irN) scattering is a fundamental hadronic process for which a large amount 
of data is available and it is important to understand this as completely as possible. Several 
relativistic models exist which provide reasonably good fits to the experimental phase 
shifts. These models consider the N, ir and A-resonance fields, the isoscalar-scalar (in some 
cases implictly via a power series expansion), the p meson and sometimes higher resonances, 
although these play a minor role. Our interest here is to examine whether a model which 
contains the minimal number of fields, namely the N, n and A , can yield equally good fits. 
Thus we effectively integrate out any other fields. For example, provided the center-of-mass 
(cm.) energy is not too high, we can expand the p propagator as (m 2 p — t)^ 1 = m~ 2 (l + 



1 



t/w? + •••), where the Mandelstam variable t = (q — q') 2 and q and q' are the initial and 
final pion cm. four- moment a. The series of terms can be absorbed into contact interactions 
in the lagrangian and it is clearly important to employ the most general set of such contact 
interactions which is consistent with the symmetries of quantum chromodynamics. 

While the A degree of freedom plays an important role in nN scattering, the Z parameter 
that specifies the form of the 7riVA vertex has been controversial, see the discussion of 
Benmerrouche et al. ||. Most of the papers cited above fit the Z parameter to the ttN 
data. This is unsatisfactory since, as we showed recently @, the scattering is independent 
of Z if the lagrangian contains the most general set of contact terms (we demonstrate this 
explicitly below). Thus results which depend on Z indicate that the contact terms have 
been implicitly constrained, whereas it is clearly preferable to employ a general lagrangian 
and allow the data itself to impose constraints. 

We would like to employ a lagrangian which explicitly embodies chiral symmetry since 
this is known to be a fundamental symmetry at low energies. Such an approach was first 
taken by Peccei || to calculate the scattering lengths and this paper represents a modern 
extension of his work to study the phase shift data. In order to systematically enumerate 
the lagrangian we can be guided by Weinberg's power counting arguments ||. For this 
purpose we identify a generic small- momentum scale Q. This is of the order of the pion 
three-momentum or the pion mass and therefore much smaller than the scale of the nucleon 
or the delta mass. Then according to the power counting, a Feynman tree diagram without 
loops contributes to ttN scattering at order Q u with 

v= l+Y i V i (d i + ±n i -2) , (1) 

i 

where Vi is the number of vertices of type i characterized by rii baryon fields and di pion 
derivatives or factors. This suggests that we associate di + \rii powers of Q to a term of 
type i in the lagrangian |ID| . Also, Krause (11]] argues that ip — M is of 0(Q), as is a single 



factor of 75 (note 7^75 is of 0(1)). Although we naively count 75 as 0(Q) for organizing 
the lagrangian, we shall show later that this counting is not precise. Chiral symmetry 



(SU(2) <g)SU(2)), Lorentz invariance, and parity constrain the possible irN interactions and 



these can be found in Ref. |fT2| . For interactions involving the A isobar we use the notation 
of our previous paper |7j and follow the discussion therein. We write the lagrangian up to 
quartic order as the sum of order Q 2 , Q 3 , and Q A parts: £ = £2 + £3 + £4 • 
The order Q 2 part of the lagrangian is 

£ 2 = N{ip + gA^l,a, - M)N + \flti {d^d»U) + \m 2 Jlti (U + tf- 2) 

+AXM + h A (A„.a v WN + NQ^a,-A u ) + K^fa** . (2) 

where the pion field arises in U(x) = exp(2m(x) / f^) with rr = |7r ■ r and the axial vector 
field a M = d^n/ f n + ■ • ■ , while the vector field = — d^n]/ f 2 + • ■ ■ . The trace is 
taken over the isospin matrices and the covariant derivative on the nucleon field is T>^N = 
d^N + iv^N. As regards the A, the kernel tensor in the kinetic-energy term is 

A"" = -(jp - M A )g^ + %{^V V + 7^) - Y{ip + M A )Y ■ (3) 

Here we have chosen the standard parameter A = — 1, because it can be modified by redef- 
inition of the A field with no physical consequences |T3|]. The covariant derivative is 



V^A V = d„A v + iv„A v -v^x A„ , (4) 

in which A M = TA M , with T a the standard 2x4 isospin | to » transition matrices. The 
off-shell Z parameter appears in Mi , = g^ v — (z + |^7 M 7i/ . We have simplified the 7rAA 
interaction in Eq. (Q) by choosing the physically irrelevant parameters Z 2 = — \ and Z 3 = 
(see Ref. @); this term does not contribute to the scattering amplitude at tree level. 
The order Q 3 part of £ is 

£3 = jpfNti {dJJ^U) - ^Nv^N 

+ ^2 m ^ V v Ntr (aV) + ^mlNN tr (U + - 2) + • • • , (5) 

where the dots represent terms that do not contribute to the nN scattering amplitude up 
to 0(Q 3 ) and we have defined 



V P = V ll - (d„ -iVf,) , (6) 
We have also applied naive dimensional analysis |14| to factor out the dimensional factors 



so that the parameters are expected to be of order unity. 
Finally, the order Q A part of £ is 

£ 4 = ^mlN^U - U^)N + ^NYD»v, u N 

+ J^ m ^{V»MT a Ntr(T a iDW,a»}) + ••• , (8) 

where the braces denote an anticommutator and 

D p a u = d p a u + i[v^ a v ] , D a v fiu = d a v fiu + i[v a , v^] . (9) 

Again the dots represent terms that do not contribute to the ttN scattering amplitude up 
to 0(Q 3 ), such terms include the usual fourth-order pion lagrangian. 

Using the pion and nucleon equations of motion |],|I^,|nj, we have simplified the contact 
terms listed in Ref. [15]. For example, we reduce the 0(Q 3 ) term iV T>^D U iVtr (a^a u ) to the 



sum of the 0(Q 3 ) K\ term, the 0(Q 4 ) A4 term, and higher-order terms which we omit. As a 
result we have the minimum number of independent terms contributing to the nN scattering 
amplitude up to 0(Q 3 ). As we have remarked, the isoscalar-scalar and isovector- vector p 
fields as given in Ref. [1(J have been integrated out. Their effects show up in the contact 



terms n 2 and A 2 . For example, in terms of the pirn coupling (g pnw ) and the pNN coupling 
(<7 P ), the rho gives a contribution to the A2 parameter of — 2g piT7 ,g p M 2 fl/rn^. 

In Fig. |I]we show the tree Feynman diagrams for ttN scattering. The crossed diagrams for 
Figs. |I](b) and ]l](c) are suppressed. The lagrangian £ 2 gives contributions to the T matrix 
of 0(Q) from all three diagrams; note that the contact diagram is due to the Weinberg 
term —N'-j^v^N. The interactions in £3 and £4 (except for the Ai term) give further 
contributions to Fig. [3](a) of order Q 2 and Q 3 , respectively. In Fig. ]l](b), each vertex can be 
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FIG. 1. Tree Feynman diagrams for irN scattering: (a) contact terms, (b) nucleon exchange, 
(c) A exchange. Crossed diagrams are not shown. 



either a pseudovector g& vertex or a symmetry-breaking Ai vertex. As mentioned earlier, 
the appearance of 75 renders the Ai vertex of higher order than expected from the chiral 
counting of Eq. flU) so we have included it in £4. The reason for this extra power of Q results 
from the following relation: 

_/a 1 / \ u{p')4u(p) 

u{p h5 tfTJ^ l5u{p) = " (p + ^-m^ ' (10) 

where u(p) is the positive-energy free Dirac spinor. Thus, with one Ai and one qa vertex, 
Fig. |l](b) is of 0(Q 3 ); we include this contribution. With both vertices of Ai type the result 



is of 0(Q 4 ), whereas, associating an extra factor of Q with each 75 as suggested by Ref. fTl 
and Eq. flD, we would expect 0(Q 5 ). 



We follow the standard notation of Holder JTTj and Ericson and Weise |18] to write the 
T matrix as 

T ba = (7r 6 |T|7r a ) = T + S ab + l[n, r a ]T- , (11) 

where the isospin symmetric and antisymmetric amplitudes are 

T ± = i4 ± + i0 + /)B ± . (12) 

Here A and are functions of the Mandelstam invariant variables s = (p + q) 2 , t, and 
u = (p — q') 2 , where p is the initial nucleon cm. momentum. They are given by the sum 



of the contributions from the contact terms in Fig. 0(a), the nucleon exchange in Fig. 0(b) 
and the A exchange in Fig. |l|(c). The amplitudes arising from the contact terms are 
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(13) 
(14) 
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where v = (s — w)/4M. The contributions from the nucleon- and A-exchange are well-known 



(see Ref. [|17| for example). We list these contributions in the following for completeness. 
The amplitudes arising from nucleon exchange are 
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where v-q = (t — 2ml) /AM. The amplitudes arising from A exchange are 
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where v A = (2M A - s - u)/4M, E A = (Ml + M 2 — m 2 )/2M A , and 
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Y(Z) 



ax = 2{E A + M)[M A (2E A - M) + M{E A - 2M)] 

M A J 



2 + mJ z+{] 



(25) 
(26) 



Notice that in agreement with Ref. J71 only the nonpole terms in the A-exchange diagram 
involve the off-shell parameter Z. Therefore these contributions can be absorbed into the 
parameters of the contact terms according to 
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We shall quote parameters obtained with Z — — | and the parameters for other values of Z 
can be obtained from Eqs. ( p7|) to (0). We have verified this numerically. 

We use the standard labelling for isospin-spin partial wave channels, namely a = 
(I, 21, 2J) where I is the orbital angular momentum, / is the total isospin, and J — I ± | is 
the total angular momentum. The elastic scattering amplitude 

1 



f a = ^e i&a sin 5 a 
\Q\ 



(32) 



is obtained from the amplitudes A ± and B ± by the usual partial wave expansion 1 19| . Here 
S a is the phase shift of the a partial wave. 

Unitarity requires f a to take the complex structure in Eq. (0). However, f a is real 
in a tree-level approximation to the scattering amplitude. We may recover unitarity by 
obtaining the phase shifts from two common methods. The first assumes that the calculated 
f a is simply the real part of Eq. (|3"2"|). The second introduces a K matrix given by |TE 



fa 



Kg 

1 — i\q\K a 



where K n 



tan 5 n 



(33) 



The calculated real tree- level amplitude f a is then assumed to actually be K a , which is true 
for \q\ small enough. For sufficiently small phase shifts, the two methods yield the same 
answer because sin5 a « tan5 a w S a . However, near the resonance region where 5 a ~ 7r/2, 
the .ff -matrix method is preferred for the following simple reason. (We note that Goudsmit 
et al. have proposed a justification for the X-matrix method.) 

First, for energies near a resonance, the amplitude in the resonant channel takes the 



relativistic Breit-Wigner form. Taking the P33 channel as an example, we have [fl£ 



Mi- s- iM A F A 
where Ta is the A width. Eqs. (|32|) and fl34 ) lead to 



Q\JP33 - „ ;u r ' ^ 



tan£ P33 = A ■ (35) 



Next, we expect that the tree- level amplitude can be obtained by setting the imaginary part 
of the denominator of Eq. (p4]) to zero: 

l«WS = ^, (36) 



and this is indeed obtained by retaining only the pole contribution of Eqs. fl2"Ip to fl2~4] ) and 
using the partial wave expansion. Finally, given the tree amplitude Eq. (^), the correct 
phase shift of Eq. ( j35|) is obtained by the f^-matrix method. Thus, while the two methods 
do not differ for small phase shifts in the nonresonant channels, the i^-matrix method is 
also good on resonance. We therefore use the ii'-matrix method here. 

In our calculations we choose the standard values M = 939 MeV, Ma = 1232 MeV, 
and m n = 139 MeV. We also take |20| f n = 92.4 MeV from charged pion decay, qa = 1-26 



from neutron (3 decay, and Ha = 1-46 from the A width, Ta = 120 MeV; allowing and 
hji to vary does not improve the fit. We first consider an 0(Q 2 ) approximation to the T 
matrix which neglects £4. The four parameters listed in Table 1 were obtained by a x 2 fit 
to the data of Arndt pD] for pion cm. kinetic energies between 10 and 150 MeV. Because 



negligible error bars are given in the data at low energies, we assign all the data points the 
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same relative weight. In Fig. |], we plot the calculated S- and P-wave phase shifts (dashed 
curves), along with the data to which we fit, as a function of the pion cm. kinetic energy; we 



also display older data from Bugg [T^ and from Koch and Pietarinen [23| . The calculation is 
in good agreement with the data up to 50 MeV, but beyond this energy the fit deteriorates 
for three of the partial waves. The value of x 2 is unity for a relative weight of 15% which is 
a measure of the accuracy of the fit. The threshold (vanishing pion kinetic energy) S-w&ve 
scattering lengths (a2j) and the P-wave scattering volumes (0,21 2 j) are given in Table f|. 
The difference between the data from Refs. |yj and |23 gives an indication of the error in 



the absence of a more reliable estimate. As regards theoretical predictions, apart from 
which is closer to the older value [23||, the 0(Q 2 ) results agree nicely with Ref. |21] which is 



to be expected since they are the zero energy extrapolation of the data we have fitted. 

We now include £4, which involves five additional parameters (Ai to A5), to take the 
tree approximation to 0(Q 3 ). The results are indicated by the solid curve in Fig. 2 which 
gives a good fit (with a relative weight of 8% for y 2 = 1) out to 150 MeV. In fact only the 
Su and P13 phase shifts deviate significantly from the data in the 150 -200 MeV range. Of 
course the rather precise agreement for Sp^ is strongly influenced by the phenomenological 
.fT-matrix unitarization. This forces the phase shift to be ir/2 at s = M\ corresponding to a 
cm. energy of 127 MeV. As regards the threshold results given in Table II, the predictions 
are a little closer to the data than at 0(Q 2 ) with the exception of a±. In this connection 
it is instructive to examine the isoscalar and isovector S wave scattering lengths, (60, &i)- 
A recent determination |24| gave (—0.008 ± 0.007, —0.096 ± 0.007) in units of m' 1 , in 



substantial agreement with Refs. |^T],^3[; note that Arndt favors a value of bo consistent 
with zero. At 0{Q 2 ) we obtain (0.007, -0.081) and at 0(Q 3 ) (-0.010, -0.077). Thus the 
isoscalar b , which is zero in the chiral limit, has improved by going to 0(Q 3 ), while the 
magnitude of b\ remains too small. 

Apart from A3 which has little influence on the fit, the 0(Q 3 ) parameters listed in Table | 
are of order unity although Eqs. fl2"?D to ( |5T] ) show that, while the fit is independent of Z, 
the actual parameter values will depend on Z. The pseudoscalar coupling with parameter 
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TABLE I. Parameters from fits to the S- and P-wave phase shifts. 



fit 






«l 




Ai 


A 2 


A 3 


A4 A5 


0(Q 2 ) 


-0.1960 


0.5001 


0.3061 


-0.9328 










0(Q 3 ) 


-0.1376 


0.5301 


0.7431 


-0.5799 


0.3650 


-0.3239 


-0.0401 


0.6334 -0.4347 



Ai allows the effective ttNN coupling constant to be adjusted in the 0(Q 3 ) fit. From 
the Goldberger-Treiman relation, our values for g& anf f n correspond to a ttNN coupling, 
9ttNN =12.8 which is a little lower than the value of 13.1 obtained by Arndt et al. |[25|| . When 
the Ai term is included g^NN decreases slightly to 12.6. We will not comment on the sigma 
term since this requires extrapolation to the unphysical region which may not be reliable 
with this tree-level model. 

With 9 parameters our 0(Q 3 ) calculation deviates from the data only beyond 150 MeV 
cm. energy. At the higher energies we do a little better than Goudsmit et ai who have 
7 parameters and fit to 75 MeV. The calculation of Boffinger and Woolcock , which is an 
improved version of Ref. contains 10 parameters and produces a fit which is similar to 
ours but a little better at energies ~ 200 MeV. The remaining models HQ have a larger 
number of parameters (14) and correspondingly fit to significiantly higher energies. 

In conclusion, we have discussed a chiral lagrangian involving just the basic N, n and A 
fields, with a series of terms representing a momentum expansion. We find that a tree-level 
calculation with this model represents the data as well as other models with a similar number 
of parameters. Further we have confirmed by explicit calculation that the Z parameter of 
the 7rNA vertex is irrelevant if a sufficiently general lagrangian is employed. Of course it 
would be more satisfactory if a unitary scattering amplitude emerged naturally, rather than 
being imposed phenomenologically. Such would be the case if loops were calculated in heavy 
baryon chiral perturbation theory and work in this direction is in progress. 

We acknowledge support from the Department of Energy under grant No. DE-FG02- 
87ER40328. 
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FIG. 2. The calculated S- and P-wave phase shifts as functions of the pion cm. kinetic energy. 
The phase-shift data from Arndt |2l| (triangles), Bugg (2^] (squares), and Koch and Pietarinen 
p3| (circles) are also shown. 
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TABLE II. The calculated /S-wave scattering lengths and P-wave scattering volumes for the 
0{Q 2 ) and 0(Q 3 ) fits compared with the data of Refs. [21| and 23]. The scattering lengths and 
volumes are in units of m" 1 and m~ 3 respectively. 



length/ volume 



0(Q 2 



0(Q 3 



Ref. [21 



Ref. 123 



ai 

03 
an 

ai3 

031 
0133 



0.169 
-0.074 
-0.074 
-0.032 
-0.038 

0.212 



0.144 
-0.087 
-0.071 
-0.031 
-0.040 

0.209 



0.175 
-0.087 
-0.068 
-0.022 
-0.039 

0.209 



0.173 
-0.101 
-0.081 
-0.030 
-0.045 

0.214 
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